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ABSTRACT 


U. S. Naval operations in the Arctic require an effective way to predict the move¬ 
ment and behavior of sea ice. This is currently provided by the Navy's PIPS model 
which is based on Hibler's (1979) sea ice model which combines the thermodynamic ice 
heat budget with a dynamic ice model sensitive to the effects of ice thickness and ice 
strength. The PIPS model simultaneously solves a system of four equations of ice mo¬ 
mentum balance, ice rheology, ice thickness, and ice strength. In order to test the per¬ 
formance of Hibler's formulation, another version of it, developed by Lemke et al., 
(1990), was adapted to the Arctic Ocean. The model was initialized and run using 1986 
forcing data and its performance evaluated using Arctic buoy drift data. Results indicate 
that the model ice drift is principally driven by wind forcing, that its response to changes 
in weather is rapid and essentially correct, and that it performs better at high wind 
speeds than at low wind speeds. Limitations to its accuracy were chiefly the result of 
limits to the precision and resolution of the input data provided to run the model, espe¬ 
cially near the ice margins. Overall, the model performs well in depicting the ice flow’ 
pattern in all conditions in the Arctic. 
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I. INTRODUCTION 


Ever since man has sailed the cold waters of the North Atlantic and the Arctic he 
has had a desire and a need to determine the extent of the ice cover, ice motion, and ice 
thickness. Modem day operations in the Arctic, especially by the U. S. Navy, have made 
it important to be able forecast the behavior of the ice cover, namely its concentration, 
thickness, and motion. Real time knowledge of the sea ice extent and ice movement are 
provided on a limited basis today by satellite imagery and sparse buoy data. No tech¬ 
nique is yet available to routinely provide real time ice thickness measurements. Thus 
polar scientists must rely on ice forecast models to provide knowledge of the temporal 
and spatial variability of these parameters. The development of effective models which 
can predict the ice concentration, thickness, and velocity is important for the long term 
study and safe use of the Arctic as well as the Antarctic. This has led to the development 
of a standard Navy model, the Polar Ice Prediction System (PIPS), which was designed 
to accomplish these goals. 

In the Arctic considerable interaction between the atmosphere, oceans, and the ice 
layer takes place requiring the dynamic and thermodynamic relationships to be quite 
complex. Naturally expressions for these relationships have evolved over time, becom¬ 
ing more sophisticated as more measurements of various ice parameters has shed more 
light on their physical behavior and interaction. Early models used empirical relation¬ 
ships to determine ice motion or drift. The first model used by the Navy was developed 
by Skiles (1968). It related ice drift to the geostrophic wind and mean upper ocean 
currents. This model was subsequently replaced by the free drift ice model developed 
by Thorndike and Colony (1982). This model was based on a relationship between the 
geostrophic wind, ice, and ocean currents determined from a statistical analysis of five 
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years of drifting buoy data. These two models had serious deficiencies in that they did 
not take into account the effects of ice thickness, ice cc ncentration, ice growth, and 
internal ice stress on ice drift. 

At the same time these emperical ice drift models were in use, sophisticated sea ice 
models were developed that concentrated on the thermodynamic treatment of ice and its 
interaction with the atmosphere and the ocean (Maykut and Untersteiner, 1971; 
Semtner, 1976; Parkinson and Washington, 1979; Manabe et ah, 1979). Early work re¬ 
lating ice dynamics and thermodynamics was conducted by Thorndike et al. (1975), 
Rothrock (1975), and Coon (1974). The concepts of an ice thickness distribution func¬ 
tion and the coupling of ice thickness to ice rheology in the studies of ridging arose from 
these early studies. In 1979 these concepts of ice dynamics were combined with a 
working thermodynamic ice model to produce the Arctic Ice Dynamic Joint Experiment 
ice model (Pritchard, 1980). Hibler (1979) developed a model employing these concepts 
along with an Eulerian coordinate system and a viscous-plastic ice rheology which made 
long term simulation feasible. In time, a version of this dvnamic-thermodyamic sea ice 
model was evaluated and modified to run at the Fleet Numerical Oceanography Center 
(FNOC) and was designated the Navy's standard model, known as the Polar Ice Pre¬ 
diction System (PIPS), for predicting ice motion, ice concentration, and ice thickness 
(Preller, 1985). 

Hibler's model, the basis of PIPS, has been in wide use since 1979. Various studies 
have sought to define the model's accuracy by comparing different aspects of the model 
results to real world ice conditions. These include studies by Hibler and Walsh (1982) 
on seasonal variations, and a study by Tucker (1983) which compared the model ice drift 
results to buoy data in the Greenland Sea. Later a study was conducted by Tucker and 
Hibler (1986) which compared the PIPS forecr u accuracy with a free ice drift model and 
buoy data that resulted in several improvements to the operational model (Preller and 
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Posey, 1989). Owens and Lenie (1990) developed a version of Hibler's ice model cou¬ 
pled to Lemke's one dimensional mixed layer ocean model for studies of the Antarctic. 

The purpose of this study was to adapt the Owens-Lemke version of the Hibler sea 
ice model and run it on the NPS mainframe computer using historical input forcing data 
and to examine its ability to describe ice motion and ice concentration as compared to 
buoy trajectories and ice edge reports in the Arctic. The forcing data came from the 
FNOC Naval Operational Global Atmospheric Prediction System (NOGAPS), a model 
that provides analysis and prediction fields of atmospheric forcing variables from the 
year 1986. Buoy trajectories and satellite derived ice edge reports from 1986 were used 
for comparison of model results. The study concentrates on the ice drift velocities in the 
central Arctic basin and ice edge position in the eastern Arctic region of the Greenland 
and Barents Seas. The results are studied to see how the model handles the effects of 
seasonal variation, periods of maximum and minimum ice extent and thickness, periods 
of ice growth and decay, significant changes in the input forcing data from day to day 
such as large storms, and the effects of a snow layer. This study will also look at an 
instance of the Beaufort Gyre reversal. 

As noted earlier, studies on the PIPS model's prediction ability were previously done 
by Tucker and Hibler (1986). During their study the model was updated on a weekly 
basis. Updating is not employed in this study. Instead the ice model's ability to portray 
conditions as accurately as possible, i.e., the real ice conditions over the course of a year, 
will be studied with historical analysis fields. This in effect removes one source of error 
(forecasting atmospheric forcing fields) from the study. It is expected that this will re¬ 
veal some strengths and weaknesses in the ice model. 

The following chapters contain an overview of the theories behind the dynamic 
thermodynamic sea ice model, followed by some specifics on the model version used in 
this study. This is followed by a chapter which discusses the input data and parameters 
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used in the ice model. The fourth chapter provides and discusses the results of the study. 
It is followed by a summary and conclusions chapter. 
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II. MODEL DESCRIPTION 


A. GENERAL DESCRIPTION 

The ice model used in this study is essentially the dynamic-thermodynamic sea ice 
model developed by Hibler (1979). As previously mentioned in Chapter I, this model 
offers improved dynamics over previous sea ice drift empirical models (Skiles, 1968; 
Thorndike and Colony, 1982) and added an ice thermodynamics component. Succinctly 
stated, the dynamic-thermodynamic ice model is an interaction between ice momentum 
balance and the continuity equations of ice thickness and compactness (Stossel et al., 
1990). A detailed discussion of the model is provided by Hibler (1979) and in subsequent 
articles (Hibler, 1984; Preller, 1985; Lemke et al., 1990). An overview of the model, 
primarily derived from Hibler's 1979 article, is provided here. The main components of 
the model include equations to determine the momentum balance, the ice rheology, the 
ice thickness distribution, and the ice strength of an ice field in a specified domain. 

Ice momentum balance describes the ice drift through wind and water stresses, 
Coriolis force, internal ice stress, inertial forces, and ocean tilt. This momentum balance 
can be discribed by the following equation 

m = — mfk xu + tj + T,- mgVH 4- F (1) 

(1) (2) (3) (4) (5) 

where the left hand side of the equation is the ice acceleration due to temporal and 
spatial changes in ice velocity (u). The right hand side of the equation represents the 
sum of the forces acting on the ice field; (1) is the Coriolis acceleration, (2) is the forcing 
due to wind stress, (3) is the forcing due to water stress, (4) is the forcing due to the tilt 
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of the sea surface, and (5) is the forcing from internal ice stress. In equation (1) m is the 
ice mass per unit area, f is the Coriolis paramtier, g is the acceleration due tdigravity, 
and H is the sea surface dynamic height. Figure 2.1 schematically shows how the forces 
act on an ice floe and their relative magnitude. 

Observations have shown that the force balance is dominated by the air and water 
stresses, the Coriolis term, and, under compact ice conditions, the ice interaction or 
internal stress term (Hibler, 1984). Measurements of water and air stresses show them 
to be of the order of 0.1 Nrrr 2 . Wind stress is the larger of the two stresses in the central 
polar pack and plays the largest role in the movement of the ice. The Coriolis force is 
approximately 0.05 Nnr 2 -, the internal ice stress is variable but can reach similar magni¬ 
tudes (Hibler, 1984). 

The air and water stresses are defined by the following equations 

T a = Pa C a I Vg I ( U g COS + k X U g sin 0) (2) 

T w = p w C w | U w - u | [((7 W - u) cos 6 + k x ( U w - u ) sin 0] (3) 

where U t is the geostrophic wind, U. is the geostrophic ocean current, C„ and C. are the 
air and water drag coefficients, p, and p„ the air and water densities, and 0 and 6 are 
boundary layer turning angles for air and water, respectively. A turning angle is the 
angular difference between the ice motion and the direction of the force acting upon it 
(in this case geostrophic winds and oceanic currents), the difference arising from friction. 
The geostrophic winds act above the thin boundary layer that overlies the ice surface. 
Water stress acts in a manner opposite to the direction of ice movement (McPhee, 1979). 
Figure 2.2 illustrates the turning angle principle. 

Ice is not completely solid; it can deform and flow like a fluid under stress. Ice 
rheology is the science of ice deformation, crushing, and flow through elastic, viscous, 
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and plastic means. Plastic behavior defines a permanant change in shape or deformation 
as a result of stress but without breaking. Viscosity determines a fluid's resistence to 
flow, while elasticity is an object's ability to return to its original shape after undergoing 
a stress. Ice breaks apart easily under tension resulting from diverging forces. During 
periods of convergence it may fracture or crush, rather easily when large amounts of thin 
ice are present. Thick ice, on the other hand, exhibits considerable resistence to com¬ 
pression and shearing. Hibler defmed the deformation stresses involved through a 
viscous-plastic constitutive law which relates the ice stress to the strain rate and ice 
thickness (or equivilantly, ice strength). The force resulting from the internal ice stress 
is defined as follows: 


F=Va (4) 

where a, the two dimensional stress tensor, is calculated by the following equation 

a (l = 2»/ (e tJ , P)e (J + [(C(e</, P) - P))c kk - — ]<5 iy (5) 

where t lt is the strain rate tensor, P is a pressure term related to ice thickness, e w is the 
ice divergence and £ and r\ are the non-linear bulk and shear viscosities, respectively. In 
this model sea ice is considered to be a non-linear viscous compressible fluid that inter¬ 
acts in a rigid-plastic manner. This allows the ice to resist compression and shear de¬ 
formation at high stresses while allowing dilation to occur when little or no stress is 
applied. In a manner proportional to the ice thickness or strength, ice is considered to 
flow plastically for high negative strain rates (compaction) and to deform in a linear 
viscous manner for sr. ill negative strain rates. This allows for ridging under com¬ 
pression and thinning under divergence. Figure 2.3 depicts how the ice behaves under 
stress. 
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The thickness of ice affects how the ice drifts and how it will deform under a stress. 
Thin ice with its relatively low mass will tend to be moved more by surface winds and 
other forces and will fracture and crush easily when compacted. Thick ice, with more 
mass per unit area, is more resistant to movement from forcing. Thick ice is also able 
to withstand higher amounts of convergence than thin ice. In the model an ice thickness 
distribution function accounts for the change of ice thickness and concentration due to 
growth, ablation, advection and deformation of sea ice. The original Hibler (1979) ice 
model considered the ice within a given grid cell to have only two thicknesses: thick ice 
of thickness h and having a concentration or compactness A (0 <, A <. 1) and thin ice of 
thickness essentially zero (includes open water) and a concentration of (1 - A). Later 
models, including the version used in this study, updated this two level thickness 
parameterization to a multilevel thickness parameterization (Hibler, 1984) with several 
levels of thickness h equally spaced from zero to twice the average ice thickness. The 
model used in this study uses a seven level parameterization. 

Growth rates for thick ice (> 0.5 m) are considerably slower than for thin ice (< 0.5 
m) (Owens and Lemke, 1990). This is a result of the slow transfer of heat through ice 
over a non-linear temperature gradient from the warmer ocean below to the cooler at¬ 
mosphere above the ice layer. Thin ice has an almost linear transfer of heat across its 
layer allowing for faster growth in the presence of freezing conditions. Since the mass 
for thin ice is small, the average thickness for both thick and thin ice in a model cell is 
considered to be h/A. The seven level ice thickness, h, and its concentration, A, are de¬ 
rived from the continuity equations of thickness and compactness and include the effects 
of advection, convergence, divergence, thermodynamic growth and decay, and diffusion. 

* A „ (wA) „ (vA) . „ . JtjBP . 

8 ~ = ~ d —— 5 — + S * + fusion ( 6 ) 
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(7) 



„ (uA) „ {vA) 

8— - 8^r- 

ox oy 


+ S A + diffusion 


where the left hand side of equations (6) and (7) are the temporal changes of ice thick¬ 
ness and concentration, respectively. The first two terms of the right side of the 
equations describe the effects of advection, and convergence/divergence with u repres¬ 
enting motion in the x direction and v represeting motion in the y direction. With 
advection ice can be moved into and out of a grid cell resulting in convergence and di¬ 
vergence, respectively. Divergence will decrease the ice thickness in a cell, opening leads 
and polynyas in which thin ice forms. Convergence will close leads and polynyas, 
crushing thinner and weaker ice thereby increasing the ice thickness in a grid cell. The 
third terms, S* and S A , are the thermodynamic terms of growth and decay. They are 
described in the following equations 

S„-A-j)A + (l- 1 4)A0) ( 8 ) 

(f-T-O-'O if AO) >0) r o ./s.xn 

5 -H ° M A f ( 9 > 

l o if m < oj >f s h < oj 


where ffh) is the ice growth rate for each level of ice thickness h and ffO) is the growth 
rate of thin ice, h„ is a fixed demarcation between thick and thin ice of 0.5 m, and (1-A) 
is the fraction of open water. The S k term is the sum of the ice grown in both open water 
and the additional growth of thick ice. 

The growth rate of ice, ffh), is computed from a complex heat budget which deter¬ 
mines the oceanic and atmospheric heat fluxes through the ice which arise from the 
conservation of energy at the upper and lower surfaces of the ice. Most ice models (in¬ 
cluding the one in this study) include provisions for calculating a snow layer for use in 
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the heat budget. The snow layer has the conflicting effects of providing an insulating 
cover to the ice surfaces which reduces ice growth rates and also having a high albedo 
which reflects incoming solar radiation reducing the rate of ice melt. The heat fluxes 
involved include incoming solar (short wave) radiation, net long wave radiation, sensible 
heat, latent heat, and conduction of heat from the ocean to the underside of the ice, from 
the upper ice surface to the snow, and from the snow layer to the atmosphere. Figure 
2.4 provides a profile of the heat budget across the ice and snow layers. Following 
Semtner's (1976) formulation, heat is transferred through the ice sheet by assuming a 
linear temperature profile within the ice and a constant ice conductivity. When heat is 
lost from open water to the atmosphere, the ice growth rates resulting from the heat loss 
are taken be vertical in direction. When open w'ater absorbs heat from the atmosphere, 
the heat is allowed to mix under the flows to reduce growth rates. Any remaining heat 
can be used to raise the mixed layer temperature, which causes lateral melting, or be used 
to raise the temperature of ice. In the presence of an ice cover the mixed layer is held 
at the freezing temperature until all the ice in a grid cell is melted first. (Preller and 
Posey,1989). 

The strength of a substance is a measure of its ability to resist deformation as a re¬ 
sult of stress. In this model ice strength (P) in equation (5) is treated as a function of 
the ice thickness distribution and concentration and is given by the equation 

P = P*h exp [ — C(1 — AJ\ (10) 

where P* is a constant that is a measure of the maximum ice strength, the point where 
compressive stresses increase to a level sufficient to cause ice to deform as a plastic vice 
a linear viscous fluid (Figure 2.3.a). The constant C is the rate at which compactness 
(A) affects the ice strength. As explained below, this relationship makes the ice strength 
strongly dependent on the amount of thin ice present, but also allows the ice to 
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strengthen as it becomes thicker. When small amounts of thin ice are present in a cell, 
the larger quantities of thick ice remaining are subject to deformation (ridge building) 
giving the ice high strength due to its increased thickness. Conversely, for a cell con¬ 
taining a large amount of thin ice, mostly thin ice is deformed which yields a lower ice 
strength. 

The overall model is an interaction of the four major sections of the ice model 
(momentum, rheology, thickness distribution, strength). They are coupled and together 
form a complete system of equations where ice velocities depend upon the ice thickness 
and concentration through the ice strength term (and internal interaction) while the ice 
thickness and concentration depend upon ice velocity through advection and conver¬ 
gence. Figure 2.5 shows a general flow chart of the dynamic-thermodynamic ice model. 
Subsequent versions of this model essentially follow the same equations as the original. 
The model version used in this study is described in the next section. 

B. SPECIFICS ON MODEL USED FOR STUDY 

The specific version of the dynamic-thermodynamic sea ice model used in this study 
(termed the NPS model) was one used by Lemke, Owens and Hibler in their studies of 
the Weddell Sea and Southern Ocean (Lemke et al., 1990; Owens and Lemke, 1990; 
Stossel et al., 1990). It is similar to the original Hibler model. The model code was 
kindly provided by Breck Owens of Woods Hole Oceanographic Institution along with 
a set of Weddell Sea data to verify the model output after adapting it to run on the NPS 
mainframe computer. 

The NPS model uses a seven level ice thickness distribution vice the earlier two level 
level distribution. The PIPS model at FNOC also currently uses the seven level ice 
thickness distribution, a feature developed when it was found that the original two level 
approach produced an average ice thickness that was too thin compared to observations. 
The seven level approach resulted in an increase of ice thickness of approximately 0.5 
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m over the PIPS domain (Preller and Posey, 1989). In the thermodynamic part of the 
model, the ice growth rate is formulated as a highly nonlinear function of ice thickness 
while each grid cell is assumed to be covered by a distribution of ice thicknesses which 
are averaged to give the thickness h for that cell (Owens and Lemke, 1990). Ice growth 
rates were calculated from the seven-level thickness distribution function bounded be¬ 
tween zero (thin or no ice) and twice the average thickness following the procedure of 
Hibler (1984). 

The model code includes a one-dimensional mixed-layer pycnocline ocean model 
developed by Lemke (1987) that is coupled to the ice model growth subroutine in order 
to determine the vertical oceanic heat flux. Hibler's earlier (1979) model treated the 
ocean as a two layer fluid and specified the vertical oceanic heat flux as either a constant 
or a monthly varying value. The main feature of the Lemke subroutine is the prognosuc 
calculation of the mixed layer temperature, salinity and depth and profiles of the 
thermocline and halocline as a function of ice production and ice velocity (Stossel et al., 
1990; Owens and Lemke, 1990; Lemke et al. t 1990). Figure 2.6 depicts the prognostic 
mixed layer model. These calculations were based on conservation of heat and salt, 
potential energy considerations, and parameterization of entrainment heat fluxes 
(Stossel et al., 1990). The PIPS model does not use this prognostic mixed layer ocean 
subroutine but instead holds the mixed layer constant at 30 m depth and provides ocean 
currents and oceanic heat fluxes via the Hibler-Bryan (1984) coupled ice ocean model. 
This latter model couples a three dimensional ocean model to a sea ice model and is 
discussed more fully in a later section. In order to emulate PIPS as much as possible, 
the NPS model runs for this study held the mixed layer depth constant at 30 m and 
variable oceanic heat fluxes were inputed as monthly gridded data. Comparison of 
model runs using the Weddell Sea data demonstrated that the prognostic and constant 
mixed layer approaches showed similar results during both periods of maximum ice 
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concentration (winter) and during summer when the ice extent was at a minimum 
(Lemke et al., 1990). 

The model codes are modular in design allowing for increased flexibilty in adapting 
the model to different experimental conditions. Model parameters and forcing input 
data were read in by separate subroutines. Most of the calculations were also handled 
by separate subroutines and brought together by a main computational loop that di¬ 
rected the flow of the model at each time step. The model equations were run in the 
following order. The model parameters and initialization fields (either through restart 
data or basic initial data fields) were read in first. Then, at the beginning of the main 
loop, the forcing data were read in for that time step. Next, the ice pressure or strength 
was calculated, which was followed by the ice velocity integration to calculate the new 
ice velocities and viscosities through an over-relaxation scheme. The relaxation scheme 
starts out with a first guess for a value at a grid point (usually the previous time step's 
value). It then runs through an iteration procedure until a specified condition is met and 
that value is accepted. The goal is to force the routine to converge to a solution as close 
as possible to the true value. The next pan of the program deals with the changes to ice 
thickness and concentration due to horizontal advection and diffusion of the ice field. 
This was folio 'ed by the thermodynamic growth calculations for ice thickness and con¬ 
centration. It is at this point that the heat budget is calculated and used in the growth 
calculations. After the output of data and statistics the old values of the previous time 
step which had been retained for the present time step were replaced by the new values 
and a new time step was begun. The horizontal grid scheme was a staggered spatial grid, 
better known as an Arakawa B grid, with u and v ice velocities on the points of the grid 
cell and thickness and concentration offset to the center of the grid cell (Mesinger and 
Arakawa, 1974; Hibler, 1979). The time stepping was accomplished by a leapfrog- 
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trapeziodal scheme (Owens, 1990). Figure 2.7 depicts a - odel grid cell with the variables 
arranged to a "o" grid scheme (Hibler, 1979). 

The flexible model codes allowed for adaptation of the model to Arctic data and to 
the PIPS model domain with little change to the model code. Once the model was 
adapted to the NPS mainframe computer with the Weddell Sea data and satisfactory 
results were obtained (from comparison to runs of the same data provided by Owens), 
it was adapted for Arctic use. Parameters used in the model were reviewed to ensure 
they were appropriate for the Arctic. The PIPS grid was incorporated and the time step 
was reduced to 12 hours to match the Arctic input data sets. Most changes to the model 
occurred in the data input section. The model was run for three years to initialize the 
output using 1986 forcing data (discussed later) similar to PIPS (Preller and Posey, 1989) 
to acheive a steady state. 

C. MODEL DOMAIN 

The domain of this model is the same as that used by PIPS at FNOC. The grid as 
described by Preller (1985) is a subsection of the FNOC northern hemisphere polar 
stereographic grid illustrated in Figure 2.8. The domain of the grid covers the central 
Arctic basin, and the Barents, Greenland, and Norwegian Seas, essentially covering the 
Arctic basin from the Bering Strait to the north of Scotland. Figure 2.9 shows the PIPS 
grid superimposed on the Arctic Ocean with the land and ocean boundaries masked out 
(Preller, 1985). An averaged mapping factor is used to approximate equal spacing for 
the FNOC northern hemisphere polar stereographic grid in the region of the model do¬ 
main (Preller, 1985). The ice model grid has dimensions of 47 by 25 grid points with 127 
km spacing at all grid points. Figure 2.9 shows the grid laid over the Arctic region 
(Preller and Posey, 1989). The velocity grid is 46 by 24 with the same horizontal spacing. 

The PIPS boundary masks provided by NOARL for ice velocities, thickness, and 
outflow are also used in this model. The purpose of the masks is to flne the bounda- 
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ries of the model. The velocity masks were on a 46 by 24 grid with constant values at 
the outer boundary and all land equal to zero and ice and open water equal to one. 
Thickness masks were in a 47 by 25 grid with zero's for the outer boundary and land and 
one elsewhere. The outflow masks were similar to the thickness masks except that the 
outflow cells were also set to zero. The combination of ice thickness of one and outflow 
of zero represented an outflow condition. Outflow for this model was permitted only 
at the southern boundary of the Greenland and Norwegian seas (Preller, 1985). 

D. MAINFRAME USE 

The ice model was run on the Naval Postgraduate School's IBM 3033 mainframe 
computer. Though the model has been run on a Sun Computer workstation (personal 
communication, Breck Owens at Woods Hole) the NPS mainframe was selected for its 
easy accessibility and ability to store and process large data sets. 

The model was run on a batch processor due to its large size and the large amount 
of data involved, which required the adding of job control lines to the begining and end 
of the FORTRAN code. The increase in model size occurred because the Arctic data 
sets were about ten times the size of the Weddell Sea data sets and the model was run 
at 12 hour time steps vice once a day. This resulted in long computer processing times, 
so the model was run in segments during initialization and subsequent runs utilizing re¬ 
start records for continuity. At its best the model ran at 300 time steps per hour of cpu 
clock time. 
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Figure 2.1 An illustration of the force balance on sea ice movement. The ice velocity 
(vector v,) is acted upon by a strong wind stress (vector r„), a water stress acting in the 
opposite direction (vector t„), Coriolis force (fkxu ) acting to the right of the ice veloc¬ 
ity, a small surface tilting term (— mgV//), and a force resulting from internal ice stress 
(/ r „„), a term highly variable in magnitude and direction (lliblcr, 197')). 
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Figure 2.2 An illustration of turning angles. The vectors t„ and r„ represent air and 
water stresses, respectively. The vector i>, is the ice velocity. The angles 0 and »/> repre¬ 
sent the boundary turning angles for water and air, respectively (McPhcc, l')7‘)). 
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figure 2.3 A visual illustration of the viscous-plastic constitutive law. The top figure is 
a one dimensional illustration of the law where as a function of ice strength. I**, the ice 
behaves as a plastic under a negative strain rate and negative stress; as strain rate and 


stress near zero the ice behaves as a highly viscous fluid. There is no stress for positive 
strain rates (spreading). The lower figure illustrates the allowable stress state for a linear 
viscous rheology with either bulk or shear viscosity, and for a rigid plastic rheology with 
an elliptical yield curve. The stress states arc plotted as a function of the principle 
components of the stress tensor a (lliblcr 1979). 
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Figure 2.4.b A schematic diagram of a 3 layer model. Temperatures in 
the layers and thicknesses of the snow and ice arc predicted on the basis 
of fluxes across internal and external boundaries (Semtner, 1976). 
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Figure 2.5 A flow chart for the dynamic-thermodynamic ice model. At the top winds and ocean 
currents arc fed into the momentum equation to determine ice velocities. Thermodynamic 
fluxes arc inputed into a heat budget that influences ice thickness along with ice advcction. This 
is sent to an estimation of ice strength which is used to determine the internal ice stress which 
in turn is used in the momentum balance. In this figure C represents Coriolis force, r„ the wind 
stress, t w the water stress, F the internal ice stress, (1 ocean currents, and T the surface tilting 
term, h the ice thickness, and f^h) the growth rate (Iliblcr, 1979). 
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Figure 2.6 A profile ofl.cmkc's (1987) one-dimensional ocean mixed layer model where 
//, is the ice thickness, h the mixed layer depth, h h is the bottom layer depth, S and S s arc 
the respective layer salinities, T and 7* are the respective temperatures, and <7, and <1, arc 
the changes in salinity and temperature, respectively (Lemke ct a!., 1990). 
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Figure 2.7 The spatial arrangement of variables used in the ice model. The variables u 
and v arc the components of ice velocity, m is the ice mass, h is ice thickness, A is ice 
concentration, I* is ice strength, r. is the ice strain rate, ( and rj arc the bulk and shear 
viscosities, and I) is the spatial difference (lliblcr, 1979). 





Figure 2.8 An illustration of the standard FNOC 63 \ < 3 northern hemispheric stci cographic 
grid. The solid line boundaries show the regions used for the Arctic and IMPS domain (Prcllcr, 
1985). 
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ch is a 47 by 25 grid with a spatial resolution of 127 km. The land 
boundaries in the model (Prcller. 1985). 
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III. FORCING DATA 


A. DESCRIPTION OF THE DATA 

Four categories of forcing data are required to run the ice model. They are wind 
stress, geostrophic ocean currents, atmospheric thermodynamic variables, and oceanic 
heat fluxes. The forcing data used in this study were kindly provided by Ruth Preller 
and Pam Posey (personal correspondance, 1990) of the Naval Ocean and Atmospheric 
Research Labratory (NOARL). The data covered the year 1986 and were in a gridded 
format when received. 

The source of all the atmospheric data used in the model originated from 1986 
analysis fields of the Naval Operational Global Atmospheric Prediction System 
(NOGAPS) model in use at FNOC. NOGAPS is a global atmospheric spectral model 
which, at the time this data set was created, was a 47 wave, 18 layer resolving model 
covering the globe (current version 3.2 is an 80 wave model). Spectral models resolve 
waves vice distinct grid points in their computations via Fourier transforms. The older 
version of NOGAPS had a coarser resolution than the current model (4.5 degrees of 
longitude) and orographic details (mountains) were not well resolved. Data sparse re¬ 
gions (like the Arctic) are filled in through an optimal interpolation process. The at¬ 
mospheric model does make use of satellite imagery and buoy data (Wash, personal 
comunication). The current NOGAPS model is run at six hour intervals. At the time 
of the 1986 data set, however, it was run at twelve hour intervals. Originally the PIPS 
model was run once a day. An improvement in the ice model was noted when the model 
was upgraded to a six hour time interval to take advantage of the higher accuracy in the 
atmospheric fluxes from NOGAPS (Preller and Posey, 1989). 
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The forcing data included data necessary to derive the geostrophic winds and 
thermodynamic variables. The raw data fields included surface pressure fields in milli¬ 
bars (mb) for deriving geostrophic winds, surface vapor pressure in mb for deriving 
specific humidity, surface air temperature in degrees Celcius, solar radiation in calories 
per square centimeter per hour, and total heat flux and sensible heat flux also in the 
same units. The last three fields were used in determining the net long wave radiation. 

Winds have a large influence on ice movement. The PIPS model originally ran with 
planetary boundary layer winds (marine surface winds); however, comparison tests be¬ 
tween PIPS , free drift models (which use geostrophic winds), and buoy data showed that 
the boundary layer winds yielded less accurate results (Hibler and Tucker, 1986; Preller 
and Posey, 1989). As a result, the PIPS model was modified to use geostrophic winds 
derived from NOGAPS surface pressure fields. Geostrophic winds are used in the NPS 
model. 

Monthly mean geostrophic ocean currents (cm'sec) and oceanic heat flux es (meters 
of ice melted per month), also in gridded format, were derived from the H er-Bryan 
(1984, 1987) coupled ice ocean model. This model is a climatological model which uses 
a three dimensional baroclinic ocean model (Bryan, 1969) coupled to a two layer 
dynamic-thermodynamic sea ice model (Hibler, 1979). The sea ice model supplies heat 
flux, salt flux, and momentum exchange for the top of the ocean. The ocean model 
supplies current velocity and heat exchange information to the ice model. Mean annual 
observed oceanic temperature and salinity data were used to weakly force the ocean 
model. This allowed the ocean model to be relaxed to available climatological ocean 
data, while allowing significant adjustment in the upper ocean layer as a result of effects 
of ice-ocean interaction. Currents due to surface stress were a part of the model pred¬ 
ictions. 
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The resulting monthly ocean currents and ocean heat fluxes have significantly im¬ 
proved the ice model's ability to predict ice edge locations compared to constant beat 
flux ice models especially in the Greenland and Barents Seas (Hibler and Bryan, 1984; 
Preller and Posey, 1989). Figure 3.1 compares ice edges between a constant heat flux ice 
model and one that uses monthly oceanic heat fluxes. Preller (personal communication) 
has indicated that a three dimensional ocean model is under development for possible 
inclusion into PIPS in the future. Though winds have a much greater effect than ocean 
currents over short time scales, ocean currents can be a factor at long time scales 
(months to years) where the currents have an almost equal effect as winds for modelling 
ice movement (Preller md Posey, 1989) and at the ice edge where strong currents can 
advect heat to melt back the ice edge in a complex fashion. 

B. READING IN THE DATA 

The following section discusses how the raw data are converted and input for model 
use. Program codes for reading the data from tape and deriving the necessary forcing 
fields were provided by Ruth Preller and Pam Posey along with the raw data fields pre¬ 
viously described. Using these programs separately from the model, the data were read 
from tape and converted to a form appropriate for use in the model and sent to mass 
storage files that could be called up by the batch processing component of the 
mainframe computer. 

Surface pressure fields were converted into u and v components of geostrophic winds 
(as well as combined wind magnitude for heat budget calculations) through pressure 
gradients. A one degree Celcius correction was added to the surface air temperature 
fields to correct for an underestimate in the NOGAPS forcing fields (Posey, 1990). 
Specific humidity was calculated from the surface pressure and surface vapor pressure 
fields. Net long wave radiation was derived from solar radiation, total heat flux and 
sensible heat flux fields. Surface air temperature, geostrophic wind magnitude, specific 
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humidity, incoming solar radiation, and net long wave radiation were grouped into the 
atmospheric thermodynamic forcing terms for use in the model's heat budget routine. 
The u and v components of geostrophic wind and ocean currents were read in for use in 
the momentum balance portion of the ice model. Oceanic heat fluxes were read in to 
the ocean mixed layer portion of the model with a fixed mixed layer depth of 30 m cou¬ 
pled to the ice growth section of the model and brought into play after the atmospheric 
heat budget was determined. 

C. UNITS AND PARAMETERS 

The ice model was designed to run with all data fields in mks units. Where possible 
labels describing the units were placed in the model for better understanding for follow 
on users. Much of the forcing data had to be converted to mks units for use in the 
model and such changes are included in the program codes. Pressure fields were con¬ 
verted from mb to N lm 2 . Geostrophic winds and oceanic currents are in m/sec. The 
heat fluxes were converted to W /m 2 . 

The parameters used in this model (Table 1) can have a large effect on the ability 
to accurately portray the ice floe characteristics properly. The ice strength parameters, 
P* and C, have been demonstrated to have a significant effect on ice thickness when 
varied (Owens and Lemke, 1990). Snow fall precipitation rates can also greatly affect ice 
thickness (more about this in the next chapter). Turning angles and coefficients affect 
ice motion through calculation of ice stresses. Their values were determined through 
studies by McPhee (1979) and others. The wind drag coefficient was determined through 
studies conducted on PIPS when its wind forcing fields were converted from marine 
boundary layer to geostrophic winds. 
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Table 1. ICE MODEL PARAMETERS. 


Parameter Name 

Value 

x-distance 

127 km 

y-distance 

127 km 

time step (dt) 

4.32E + 04 sec 

Coriolis parameter 

1.46E-04 sec -1 

Density of ice p in 

0.91E + 3 kgm-' 

Density of water p w 

1.0E + 3 kg nr 1 

Density of air p„ 

1.3 kgm - 3 

Density of snow p !now 

0.33E + 3 kg nr 3 

Thick ice cutoff h 0 

0.5 m 

Ice strength P* 

2.75E + 04 Nm~ 3 

Compactness parameter C 

20.0 

Wind drag coefficient 

0.8E-3 

Water drag coefficient 

5.5E-3 

Air turning angle <t> 

25 

Water turning angle d 

25 

Reciprocal of heat of fusion 

3.31E-9 

Bulk sensible heat transfer coefficient 

2.284 

Water latent heat transfer coefficient 

5.6875E + 3 

Ice latent heat transfer coefficient 

6.4474E + 3 

Stefan-Boltzman constant 

5.5E-8 cal nr 1 °K~* 

Conductivity of ice 

2.1656 Wm~' “AT 1 

Albedo of ice 

0.75 

Albedo of water 

0.10 

Freezing temperature in the ocean 

-1.96 °C 
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Figure 3.1 Mode] simulated thickness using (a) a constant oceanic heat llux o( 2 
W -m’, (b) monthly mean Iliblcr-Bryan oceanic heat flux. I he southernmost 
contour approximates the model ice edge. Contour interval is 0.5 rn (Prcllcr and 
Posey, 1089) 





IV. RESULTS AND DISCUSSION 


The ice model was initialized with a uniform ice thickness of 3 m, an ice concen¬ 
tration of 100 percent, ice velocity of nearly zero, and air temperature at freezing over 
most of the model domain except for the southern Greenland and Norwegian Seas where 
ice concentration is zero. The model was allowed to "spin up" to a stable condition for 
a period of 3 years using the 1986 data input forcing fields described previously. At the 
end of each year period the model output was placed in a restart file for use in the next 
year's run. Analysis of the total ice volume and ice extent of preliminary runs showed 
that the model stabilized fairly rapidly in the latter part of the first year and throughout 
the second year. The remainder of this chapter is focused on the analysis and discussion 
of the results from the last (third year) run of the ice model. This involved examining 
the ice model output ice drift rates, ice extent, and ice edge concentration and comparing 
them to ground truth data for the same time period. Ice thickness, also an ice model 
output, has no independent equivalent measure for comparison so only a qualitative 
description will be given where necessary. 

A. ICE EXTENT 

The actual sea ice area within a given basin as defined by Comiso (1990) is the sum 
of the products of the area of each data element or grid cell and the corresponding ice 
concentration (in which 0 < A < 1) of that grid cell. Ice extent is defined as any area 
that contains 15 percent or greater ice concentration (Parkinson et al., 1987). Ice extent 
includes areas of open water such as leads and polynyas whereas ice area considers only 
that part which is ice. Articles by Parkinson and Cavalieri (1989), Zwalley et al. (1987) 
Comiso (1990), and Gloersen and Campbell (1988) have noted an annual oscillation in 
ice extent that is fairly stable (see Figure 4.1). During the summer the ice area is on 
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average observed to be about 78 percent of the ice extent in the Arctic basin (Comiso, 
1990). 

Ice model runs conducted during this study showed a seasonally varying ice extent 
that had a maximum of 9.41 x 10 6 km 2 at day 130 (10 May 1986) and a minimum of 
7.3 x 10‘/cm 2 on day 260 (17 September 1986) or a 22 percent change in ice extent. The 
modeled ice extent is plotted in Figure 4.2. These values are compared to those deter¬ 
mined from a study conducted by Parkinson and Cavalieri (1989) using satellite passive 
microwave energy which estimated the sea ice extent both regionally and over the entire 
Arctic region. A combination of regions including the Arctic basin, Kara, Barents, and 
Greenland Seas approximated the NPS sea ice model's domain. The Parkinson and 
Cavalieri maximum ice extent was found to be 10.15 x 10*/cm 2 and minimum ice extent 
7.59 x 10 t km 1 or a 25 percent change in ice extent. The difference between maximum ice 
extent values (0.74 x 10 ‘km 1 or 8 percent ) is a result of the difference in the size of the 
areas under consideration. The southern portions of the Greenland, Kara, and Barents 
Seas are not included in the ice model domain, thereby reflecting a smaller maximum. 
The difference between minimum ice extent values (0.29 x 10 i km l or 4 percent) is con¬ 
sistent with the difference in maximum ice extent in that the ice model domain encom¬ 
passes less area than the satellite coverage. The smaller difference during the period of 
minimum ice extent is reflected in the fact that more was ice melted outside the model 
domain than was melted inside. 

Figure 4.3 is a plot of the model ice area. On day 130 (10 May) the ice area is 
9.02 x 10*/cm 2 or 4 percent less than the ice extent. On day 260 (17 September) the ice 
area is 5.74 x lO^m 2 which is 21 percent less than the ice extent, similar to the 
summertime difference reported by Comiso (1990). 

Figure 4.4 is a plot of the model ice volume which is the summation of the product 
of the thickness of each grid cell and the ice area of that cell. The maximum ice volume 
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occurred on day 130 (10 May) with a value of 32.09 x 10 i km i which equates to an aver¬ 
age ice thickness of 3.55 m. The minimum ice volume occurred on day 260 (17 Sep¬ 
tember) with a value of 19.47 x 10 l km* which equates to an average ice thickness of of 
3.46 m. Discussions with Ruth Preller (personal communication) indicated that the NPS 
ice model's thickness was a little larger (up to 0.5 m of thickness difference) than similar 
runs with PIPS. A possible reason was that the NPS ice model was run without a snow 
layer which allowed thicker ice to grow due to a lack of insulation, and which subse¬ 
quently was not melted as fast as thinner, snow-covered ice might have been. Also in 
support of the snow argument is the finding by Owens and Lemke (1990), w'ho noted in 
Weddell Sea runs of the ice model, that the presence of snow substantially lowered the 
summer minimum ice extent. 

B. COMPARISON OF MODELLED ICE DRIFT TO BOUY DATA. 

The ice model drift rates or ice velocities w'ere studied in a qualitative manner under 
several different environmental conditions such as calm or light winds, strong storms, 
large broad pressure gradients, gyre reversals, and minimum and maximum ice condi¬ 
tions at selected geographic regions. The focus was to see how the ice motion behaved 
under each of these conditions. Ultimately seven cases were selected for comparison. 
Each is described below. The modeled ice drift rates w r ere compared to Arctic buoy drift 
information. These buoys, part of the Arctic Ocean Buoy Program first established in 
1979, are dropped in the Arctic ice pack, and stay with the ice until it melts or departs 
the Arctic region. Buoy postions, as well as pressure and temperature, are relayed via 
System ARGOS to a processing center w r here drift rates and pressure fields are derived 
(Thorndike and Colony, 1980). The results of this data processing are described in an¬ 
nual reports; the 1986 results are in Colony and Rigor (1989). 

Since the buoys are not uniformly distributed throughout the Arctic basin, compar¬ 
isons between the ice model and buoy drift rates were restricted to those areas of the 
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Arctic where the buoys we . present. Buoy drift rates were determined directly from the 
plots of buoy drift and surface pressure fields provided in the reports, where 1 cm was 
equal to 20 km per day or 23.1 cm/s. The ice model drift rates were measured directly 
from the vector plots of the model output data in which a 5 cm/s, 10 cm/s, and 20 cm/s 
scale was provided on each plot. General descriptions of buoy direction were made rel¬ 
ative to the surrounding geographic features and prevailing winds. 

It is important to note that a perfect correlation between buoy and model ice drift 
rates is not expected to be exact since the pressure analysis fields depicted in the buoy 
reports do not exactly match the NOGAPS surface pressure fields used to drive the ice 
model. The plotted surface pressure fields in the buoy reports were derived directly from 
the buoy measurements and additional observation stations around the periphery of the 
Arctic. The NOGAPS pressure fields are derived from a combination of data sources. 
These sources include the measurements taken from the buoys, satellite temperature 
soundings, satellite cloud winds, the previous time period's pressure analysis, and direct 
observations from ground based and air observation stations. These diverse inputs are 
weighted as to their accuracy as they are incorporated and smoothed into the NOGAPS 
analysis for the most complete and accurate analysis possible (Wash, personal commu¬ 
nication). Though NOGAPS pressure would be expected to be more accurate than buoy 
pressures alone, there might be instances where buoy-derived pressure fields pick up a 
weather feature that was missed or smoothed out of the NOGAPS processing. Gener¬ 
ally, in most cases the buoy pressure fields and NOGAPS will show similar weather 
patterns over the Arctic. 

An exagerrated or extreme example of this difference is illustrated in Figures 4.5 and 
4.6. Figure 4.5 is a plot of the surface pressure and buoy trajectories typical of plots 
from the Arctic ouoy reports. It shows that the eastern Arctic is dominated by one low 
pressure center north of Svalbard. Figure 4.6 is a plot of the geostrophic wind trajecto- 
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ries derived from the NOGAPS pressure fields for input into the ice model. Since the 
winds are geostrophic, their direction can be equated to isobars of pressure. Wftiat can 
be seen from this figure is that there are two low pressure centers (vice one) depicted by 
the NOGAPS data; with one north of Svalbard and one north of Ellsmere Island in the 
Canadian Archipelago. The presence of these two cyclones is reflected in the model ice 
drift trajectories (Figure 4.7) which also show the two cyclonic vortices in the ice move¬ 
ment. Figures 4.6 and 4.7 illustrate the importance of accurately portraying the wind 
fields for modelling ice movement since the ice flow pattern predominately mirrors the 
wind flow pattern above it. This is especially true for regions dominated by high wind 
speeds. 

1. Case One: Decreasing High Wind Conditions 

This case study examines the ice model's ability to model ice movement under 
declining wind speed conditions; specifically to see if a time lag occurs between wind 
speed changes and subsequent ice motion. The region examined considers only one side 
of the Arctic, i. e., north of the coasts of Alaska, the Canadian Archipelago, and 
Greenland where the buoys for this case study are clustered. Figure 4.8 (a,b,c,d) shows 
the prevailing wind conditions over the Arctic region and the buoy trajectories for the 
period 18-21 January 1986. Strong winds were prevalent over the region where the 
buoys were located as a result of an interaction between a high pressure center over the 
western Arctic and a low pressure center stalled over Greenland. These high winds 
steadily decreased in intensity through 21 January as the low filled and moved eastward 
out of the area. 

The buoy drift rates reflected these high winds reaching velocities as high as 24 
cm/s north of Ellsmere Island on the 18th, but declining thereafter on each successive 
day (16 cm/s on the 19th, 9 cm/s on the 20th, 3-7 cm/s on the 21st). A large drop in 
buoy speed between the 19th and 20th occurred even though the su; ace pressure gra- 
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dient was still fairly strong at that time; however, it is noted that ice was moving toward 
the coast at the time. The buoy trajectories indicated that on the 18th to 20th the ice 
was moving pretty much as a unit but by the 21st shear zones must have formed as the 
ice north of Ellsmere had come to a stop while ice near Alaska and Greenland continued 
moving. 

The modeled ice drift output for the area where the buoys were clustered during 
18 and 19 January exhibited the same type of westward ice flow north of the Canadian 
Archipelago as the buoy trajectories. Good agreement in trajectory direction is noted 
but the model indicated higher speeds (up to 38 cm/s on the 18th, 23-38 cm/s on the 
19th). On the 20th the wind backed forcing the ice and buoys to a more southerly course 
toward the coast (Figure 4.8.c) but remained essentially unchanged in the ice model 
simulation (Figure 4.9). Model trajectories were as large as 45 degrees to the right of the 
buoy drift tracks. The ice model simulation continued to follow the wind flow pattern 
shown in Figure 4.10 with modeled ice drift speeds of 20-25 cm/s near Ellsmere Island. 
By the 21st the model ice drift had changed to a more southerly direction toward the 
coast and its speed had decreased to 10 cm's (Figure 4.11). In comparison, buoys in this 
area had already stopped moving by the 21st (Figure 4.8.d), but 12 hours later (one 
model time step) the model ice drift had also stopped moving. This would seem to in¬ 
dicate a 12 hour time lag in the response of the model to decrea. ng wind speeds, possi¬ 
bly indicating that the effects of friction and ridging in reducing momemtum are 
parameterized weaker than they should be. In summary, during high winds, the model 
speeds are slightly too fast. A delay of one day in the model to shift the trajectory to 
the south vice southwest, and a 12 hour delay in response to decreasing wind speed is 
also noted. 
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2. Case Two: Light Wind Conditions 

This case study examines the ice model's performance under the absence of any 
strong wind forcing. There were several occasions for which the winds over the Arctic 
basin were calm. Figure 4.12 (a,b,c,d) shows the pressure field and buoy trajectories for 
18-21 February 1986, a period of generally calm wind conditions resulting from an ab¬ 
sence of any strong pressure gradients. The buoys were located mostly clustered in the 
same location as in the previous case but also with two buoys located north of Bering 
Strait. Generally little buoy movement was noted during this time period except near 
the Alaskan coast where buoy speeds varied from 4 to 7 cm/s. Little directional move¬ 
ment is noted. 

During this period the ice model exhibited slow movement of the ice field with 
speeds under 5 cm/s. Figure 4.13 is the model ice drift for 19 February. It shows very 
slow ice drift rates except in the vicinity of Svalbard (the result of a low pressure system) 
and a small area where ice is being driven towards the Alaskan coast. The NOGAPS 
winds (Figure 4.14) driving the ice model show wind speeds less than 10 m/s over the 
Arctic basin except as previously mentioned. The ice model drift rates for the 19th 
compare favorably with the buoy trajectory drift rates for this date; however, consider¬ 
able difference in the direction of flow between ice model and buoy trajectories was 
noted. The buoys demonstrate a westward drift in the area north of Greenland and 
Ellsmere, while the model exhibited a more northerly direction indicating a directional 
difference of 90 degrees to the right of buoy drift. The ice trajectories appear to be 
generally following the wind pattern in this location. 

In comparing the light forcing from wind fields to the model ice drift there ap¬ 
pears to be some agreement in that light winds generated very little ice motion. There 
is considerable variability in direction however, also a reflection that under lighter wind 
stress other components of the momentum equation such as ocean currents, Coriolis 
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force, and internal ice stress exert a larger effect on the ice drift than during periods of 
high wind stress. Also at lower wind speeds the turning angle between wind and ice is 
expected to be larger which would create a large difference in direction between wind and 
ice motion. 

Calm conditions also existed on 1 March, 21 March, and 31 March 1986 as well. 
The modeled ice drift speed was in agreement with that indicated by the buoy drift re¬ 
port. As before, under these light wind conditions, the buoy and modeled ice floe tra¬ 
jectories were highly variable from day to day. Considerable differences in direction of 
drift existed between the buoy reports and the model output with ice motion varying as 
much as 45 degrees to the right or left of buoy drift direction. In summary, under light 
wind conditions ice model speeds are handled well but direction is poor with disagree¬ 
ment in excess of 45 degrees. Direction did vary to either the right or left of buoy di¬ 
rect: jn. 

3. Case Three: Gyre Reversal 

One of the better known features of the Arctic basin circulation is the 
anticyclonic rotation of the Beaufort Gyre in response to the atmospheric high pressure 
center resident over it throughout most of the year. During the summer, however, the 
gyre rotation is reversed by the action of storms or atmospheric low pressure centers that 
move into the area (McLaren et al., 1987; Serreze et al., 1989). This case study looks 
at the ice model's ability to respond and depict this gyre reversal. The 1986 buoy reports 
show several instances of gyre reversal, two of which occurred on 12-17 June (Figure 
4.15) and 26 July to 1 August(Figure 4.16), respectively. The buoys are located in an 
area extending from Bering Strait to Ellsmere Island. On 12 June the buoy trajectories 
show a clockwise rotation with a drift rate of 7 cm/s off the coast of Alaska which was 
associated with an atmospheric high in the area. By 15 June a developing low had 
moved into the area and begun turning the ice counter-clockwise near the Canadian 
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Archipelago. Peak reversal occurred on 16 June with velocities reaching 20 cm/s. This 
change in buoy drift indicates that the movement of ice appears to be strongly tied to 
the large scale wind forcing. 

For the period of 12-17 June the ice model closely followed the ice flow pattern 
depicted in the buoy drift reports. The modeled ice drift for 12 June (Figure 4.17) dem¬ 
onstrated a similar anticyclonic rotation of the gyre off Alaska to the one depicted in the 
buoy report but at speeds of 2 to 3 cm/s less than the buoy drift rate. Both the model 
and buoys showed little or no ice movement near the Canadian Archipelago. By 15 June 
a developing low had begun to turn the ice cyclonically. Again the buoy reports dis¬ 
played buoy drifts that were a little faster (2 to 3 cm/s) than the modeled ice drift; the 
same was true for the 16th. On the 17th the model ice speed increased becoming similar 
to that reported by the buoys (Figure 4.18) possibly indicating a slight lag in speed of 
2-3 cm/s. At this time the gyre reversal was fully developed with the ice off the Alaskan 
coast moving at approximately 16 cm/s toward the coast. Off the Canadian coast ice 
was moving northwestward away from the coast at approximately 4 cm/s. 

Another gyre reversal took place on 26 July to 1 August (Figure 4.16). The 
buoys were distributed over a wide area of the western Arctic basin from Siberia to just 
north of Greenland. On the 26th the ice was rotating in a clockwise manner oflf the coast 
of Alaska at 7 cm/s in response to a weak high pressure ridge in the vicinity. By 28 July 
a developing low began turning the ice counter-clockwise. This pattern persisted 
through the 3rd of August with drift rates up to 14 cm/s. The modeled ice drift on 26 
July showed similar drift speeds and an anticyclonic trajectory with ice floes moving up 
to 10 cm/s off the coast of Alaska, and varying between 5 and 7 cm/s elsewhere (Figure 
4.19). The modeled ice drift did not show a complete gyre but did show a clockwise 
turning of the ice. The model matched the buoy reports in showing the counter¬ 
clockwise turning of the ice commencing on the 28th of July in the presence of a devel- 
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oping low showing no lag in response time. The ice model speeds and directions closely 
matched those of the buoys. Speeds up to 17 cm/s were indicated off the eastern 
Siberian coast by both the buoy data and the model. The model picked up the sharp 
clockwise tum in ice drift direction north of the Bering Strait which was reflected in the 
buoy reports. Around the center of the growing cyclonic gyre, model drift speeds of 12 
cm/s were measured which was close to the 14 cm/s measured from the buoy reports. 
Like the buoy reports, the modeled cyclonic gyre persisted in the presence of the storm 
showing good agreement in speed and direction. By 1 August the counter-clockwise 
(cyclonic) gyre had come to clearly dominate the central Arctic basin (Figure 4.20). The 
ice model indicated drift rates that were 2 to 3 cm/s higher than the comparable buoy 
reports but directions of movement were similar. 

The model response to both gyre reversals was virtually immediate and matched 
the movement of the buoys. The small difference in speeds between model and buoy 
was possibly a reflection of errors in the wind forcing fields. Differences in direction 
were also small but the model ice floe pattern closely matched the buoy trajectories. 
Both the model output and the buoy report for 1986 appear to indicate that the cyclonic 
nature of the Beaufort Gyre was more of a transitory feature in the ice pack rather than 
a persistent one lasting throughout the summer. The gyre shows up more clearly when 
a storm stalls over the basin permitting more energy to be pumped into the cyclonic 
circulation. In summary, the ice appears to follow the movements and direction of the 
wind more than any other factor. 

4. Case Four. Large Strong Pressure Gradient 

This case study examines the model's performance under high wind conditions 
resulting from a strong surface pressure gradient, i. e., under conditions when the wind 
velocity is already high. A subsequent case study will examine the model resp mse to a 
growing storm. 
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Examples of strong gradients were selected that covered a large portion of the 
Arctic basin to see how the model handled large scale movements of ice. One such ex¬ 
ample was during the period of 7-12 April (Figure 4.21) when a strong storm off the 
coast of Siberia had developed in the vicinity of the Kara Sea. This created a large, 
strong gradient pattern over much of the Arctic basin, especially in the eastern portion. 
The buoys at this time were mostly on the opposite side of the basin and on the pe¬ 
riphery of the storm, but its effects on the ice movement could be discerned and com¬ 
pared to model results. Prior to April 7, buoys and model output indicated that much 
of the ice pack was motionless (except at the edges) with calm wind conditions prevalent. 
On the 7 April a sudden increase in ice movement was observed as the result of the 
storm. Between Svalbard and Greenland the buoy drift speed was measured at 16 cm/s. 
Closer to Ellsmere Island buoy drift rates attained values between 7-12 cm/s, and 2-4 
era's near Alaska. The ice flow pattern paralleled the pressure contours. Figure 4.22 
shows ice model drift trajectories for the 7 April. Model ice drift rates in the vicinity of 
the buoy near Svalbard showed speeds of 14 era's with an approximately 10 degree dif¬ 
ference in direction to trie left of the buoy direction. This difference is probably the re¬ 
sult of inadequately modeling the flow of the strong East Greenland Current, i. e., the 
model uses monthly oceanic currents. North of Ellsmere model drift speeds of 11 cm/s 
were measured with a general agreement with buoy velocity direction. Modeled ice drift 
near Alaska matched buoy drift trajectories in the same area. 

By 8 April (Figure 4.23) a wider area of the ice pack was influenced by the storm 
as it moved eastward. Both the model and buoy reports showed drift rates reaching 25 
cm/s west of Svalbard, up to 14 cm's north of Greenland, and 5 cm/s near Alaska. 
Modelled ice drift direction appeared to generally agree with the buoy trajectories 
showing a piling up of ice on the north Greenland coast and a southerly flow between 
Greenland and Svalbard. By 10 April the low had begun to fill and the pressure gradient 
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had begun to weaken. The buoys and model both output show drift rates of 9 cm/s for 
ice in the central Arctic. West of Svalbard the buoy had slowed to 14 cm/4-fend the 
model output reported 12 cm/s. Near the Canadian Archipelago the buoy reports 
showed very little movement of the thick ice; however, the model continued to indicate 
a 5 to 7 cm/s movement of ice in this area. By 12 April the pressure gradient resulting 
from the storm had virtually disappeared and calm wind conditions prevailed. Most 
buoys showed little or no movement; however, the model persisted in showing some 
movement off the Canadian coast of 3 to 4 cm/s which was consistent with the earlier 
case studies which demonstrated a 12 to 24 hour lag time in response to decreasing winds 
and ice velocities. Elsewhere the model agreed with the buoy reports. 

Two other instances of strong gradients resulting from storms occurred on 27 
September and 27 October and are reflected in the buoy reports. On the 27th of Sep¬ 
tember (Figure 4.24) a strong pressure gradient extended from the Siberian coast to 
Greenland. Buoy drift rates of 14 cm/s near the high gradient region were recorded the 
with model output (Figure 4.25) for the same area exhibiting a similar ice drift rate and 
direction. In the center of the high gradient region modeled drift rates were calculated 
up to 18 cm/s. Near the Canadian Archipelago where the pressure gradient was weak 
the buoy reports showed little or no movement. 

On the 27th of October (Figure 4.24) the surface pressure gradient extended 
from the northern Scandinavian coast to Alaska, cutting a wide swath through the cen¬ 
tral Arctic basin. Buoy drift rates were measured at 12 cm/s in the middle of the high 
gradient region; the model drift rates under these high wind speed conditions were 
overestimated by 3 cm/s (Figure 4.26). Near the Canadian Archipelago the buoys and 
model were in good agreement showing little or no movement. The ice flow direction 
from the model and the buoy data were a close match. 
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Overall in these different examples of strong pressure gradients the ice model 
performed well in matching the general flow pattern of the buoy drift trajectories. Dif¬ 
ferences in drift speed between between the model and buoys was not consistent with the 
modeled speed sometimes being higher or lower than buoys in the same area. The speed 
difference was never more than 2-3 cm/s and may be the result of the spatial resolution 
of the NOGAPS forcing fields not being able to completely match the local weather 
conditions over the data-sparse Arctic. 

5. Case Five: Large storms 

This case study examines the ice model's response to a growing storm or in¬ 
creasing wind flow pattern. Case One had noted a time lag of 12 hours under decreasing 
wind speed conditions. There were several large storms in 1986 and their impact on the 
ice pack could be seen on many occasions in both the buoy data reports and the ice 
model output. Examples of storms or low pressure centers were selected that were for 
the most part entirely contained within the Arctic basin with the idea of examining the 
cyclonic ice flow circulation generated by the storm. The pressure gradients generated 
by the storms investigated in this section were generally of smaller scale and more cir¬ 
cular than those of Case Four. Storms that occurred on 16 to 19 August, 16 December, 
and 19 June will be discussed. 

During the period of 11-15 August generally calm wind conditions prevailed 
over the Arctic with only modest ice drift rates of 5 to 7 cm/s being observed. On 16 
August a large storm developed in the western Arctic basin creating a large cyclonic ro¬ 
tation of the ice (Figure 4.27). The storm's impact on the ice pack was immediate and 
this was reflected in both the buoy report and ice model output. Buoys located 
throughout the storm region displayed a cyclonic circulation recording drift rates as high 
as 18 cm/s in the strongest part of the storm and reducing to 14 cm/s toward the outside 
of the vortex. Near the center of the low one buoy was measured with a drift speed of 
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7 cm/s. The modeled ice drift rates (Figure 4.28) showed magnitudes as high as 17 cm/s 
in the same area as the buoy located in the strongest part of the storm. Near the center 
of the storm, buoy and modeled drifts indicated speeds of 6 cm/s. On the back side of 
the storm near the Siberian coast, buoy and model drift rates indicated speeds of 12 
cm/s. The ice model drift directions compared well with buoy trajectories staying within 
10 degrees of the buoy direction and usually to the left of the buoy direction. This result 
indicates the model trajectories more closely follow the pressure contours, an indication 
of the dominant influence of wind forcing relative to other forcing terms, i.e., the oceanic 
stress and internal ice stress may have been too weak. Though there were no buoys in 
the center of the low, the model correctly showed little movement at the center with 
speeds increasing toward the outside up to 22 cm/s (in areas other than where buoys 
were located). Overall the buoy drift and modeled ice drift agreed well with each other 
in both direction and speed, and showed an almost immediate response to the presence 
of the storm. 

The storm continued to dominate the western Arctic during 17 and 18 August 
with buoy and model trajectories continuing to agree. Both indicated ice drift rates of 
upwards of 23 cm/s in the cyclonic vortex set up by the storm and less movement near 
the thicker ice of the Canadian Archipelago. By 19 August (Figure 4.29) the storm had 
begun to weaken with buoy reports showing speeds of 12 to 16 cm/s with the highest 
speeds occuring near the Alaskan coast. The model also indicated speeds of 12 to 20 
cm,'s with the highest drift rates off Alaska. The ice drift directions continued to remain 
in agreement between model and buoy except for a region near Siberia where a 45 degree 
difference to the left in heading was observed between the buoy and ice model drift. The 
buoy reports showed the buoys heading toward the Siberian coast while the model tra¬ 
jectory direction indicated that the ice was still caught up in the cyclonic circulation. 
This may be similar to the calm wind condition previously discussed where the drift di- 
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rection is highly variable at low wind speeds The speed of both reports w:« 5 cm/s; a 
fairly slow speed. This difference may have been the result of a small scale weather 
feature in a data-sparse region that was smoothed over by NOGAPS but was picked up 
by the buoy drift. With the Arctic being such a data-sparse region there will always be 
a limitation on the accuracy of NOGAPS and PIPS to accurately predict ice movement. 

On 16 December (Figure 4.30) a large storm was centered just north of Canada. 
It had a large effect over much of the Arctic basin setting up a cyclonic flow pattern in 
the Canadian basin. This storm resembles some of the large gradient features discussed 
in Case Four. But, unlike the previously discussed storm, this storm generated some 
movement of the thick ice near the Canadian Archipelago. This was detected by both 
the buoy reports and the ice model (Figure 4.31). Ice velocities of up to 17 cra/s were 
indicated by both the buoy drift and the ice model. Near the Canadian Archipelago 
both also indicated a general offshore direction of the ice drift. Where the storm pres¬ 
sure gradient was strongest the buoy reports indicated speeds up to 18 cm/s and the 
model had speeds up to 20 cm/s for ice drift. One buoy near the center of the storm 
exhibited an anomalously high drift rate of 14 cm/s, a speed inconsistent with its location 
in the buoy-derived cyclonic vortex. However, it matched the model drift speed and di¬ 
rection which suggests that the storm center was correctly placed farther west by 
N'OGAPS. In the Beaufort Sea region (near Alaska) buoys flowing back toward the 
Canadian Archipelago were measured to be moving at speeds up to 12 cm/s. However, 
the model indicated slower speeds of 8 cm/s for the same area. The speed differential 
may be the model’s attempt at parameterizing the effect of piling up of ice or a second¬ 
ary anti-cyclonic circulation working to reduce the eastward drift. 

An example of a summer storm occurred on 19 June where it was centered in 
the western portion of the Arctic basin (Figure 4.30). Generally calm wind conditions 
prevailed everywhere in the Arctic except in the immediate vicinity of the storm. This 
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is reflected in the buoy trajectories and model output. Figure 4.32 shows little ice 
movement outside the immediate storm area. Inside the storm area itself, speeds up to 
24 cm's are indicated by both the ice model and buoy drift rates with good agreement 
in direction. The high speeds generated may be the result of thinner ice at lower con¬ 
centrations during this time of the year allowing less restriction in movement. 

6. Case Six: Comparison at Maximum and Minimum Ice Extent 

This case study compares model performance and behavior of ice drift at its 
maximum and minimum ice extent. The maximum ice extent occurred near 20 May and 
the minimum ice extent occurred near 17 September. The region to be examined ex¬ 
tended from Alaska to Ellsmere Island where most of the buoys were located during 
both ice extent periods. Figure 4.33 shows the buoy and pressure fields for 20 May and 
1 7 September. Though mostly light wind conditions prevailed some significant differ¬ 
ences in the weather pattern is noted between the two dates. On 20 May a weak high 
pressure center was positioned north of Alaska. On 17 September a slight low pressure 
center was found north of the Canadian Archipelago. This resulted in buoy drift tra¬ 
jectories moving in somewhat opposite directions on these two dates. 

What is interesting to note is that the wind forcing for both dates was under 10 
cm/s for the area. Even though the wind forcing was similar on these dates the magni¬ 
tude of the buoy drifts differed. On 20 May the buoys near Ellsmere Island were sta¬ 
tionary while the buoys in the Beaufort Sea were moving westward at 6 cm/s. On 17 
September the buoys near Ellsmere Island were moving eastward at II cm/s and a buoy 
near Alaska had a drift rate of 14 cm/s. This difference in speed is most likely related 
to the reduced ice concentration in summer and relatively thinner ice at this time with 
more leads and polynyas present in the ice pack. This allows for greater movement of 
the ice than would be permitted during a period of maximum ice extent. 


48 





Figures 4.34 and 4.35 depict the ice model drift trajectories for 20 May and 17 
September, respectively. On 20 May the modeled ice drift rates were very slow 
throughout most of the basin, not attaining the velocities registered by the buoys in the 
Beaufort Sea. In the absence of strong winds the model does not move the ice much 
when the ice extent covers most of the basin. Very little can be said about the drift di¬ 
rection due to the slow rate of drift other than the direction is highly variable as the 
lesser effects of ocean currents and internal stress begin to more dominantly influence 
the ice drift. On the other hand the ice model performed well on depicting the ice drift 
on 17 September. The ice was moving a lot freer matching the buoy drift rate and di¬ 
rection with speeds of 11 cm/s eastward along the Canadian coast. The model also 
matched the 14 cm/s speed and direction of the buoy off Alaska on the 17th. Generally 
the model performed better during the period of minimum ice extent when free drift 
conditions were predominant. 

7. Case Seven: Ice Miovement in the Greenland Sea 

This case examines model performance in the Greenland Sea area. No buoys 
are available for comparison however so this discussion will primarily be an observa¬ 
tional one, looking at ice flow patterns. Figure 4.36 is the modeled ice drift for 23 Jan¬ 
uary 1986. The group of high velocity vectors in the lower left represent ice flow in the 
East Greenland Current (EGC) near the ice edge jet with speeds in excess of 30 cm/s. 
Such high speeds are nearly always present as a result of the baroclinic flow caused by 
the strong East Greenland Polar Front (EGPF). Ice drift rates along the Greenland 
coast but well removed from the EGPF reach speeds up to 10 cm/s. Enhancing these 
high speeds was a storm located in the Norwegian Sea off the east coast of Svalbard 
which generated a strong atmospheric pressure gradient parallel to the axis of the EGC. 

The modeled ice drift in the Greenland Sea can be highly variable depending on 
the prevaling weather as Figure 4.37 illustrates. Here the ice along the Greenland coast 
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has slowed almost stopping with no specific direction or drift. Away from the coast near 
the ice edge, a strong flow to the north is depicted. This is an example of the model 
basically failing to depict the ice drift in this area. The effect of the EGC is completely 
missing in this case. The ice movement appears to be tied more to the forcing of the 
wind than to the ocean current. The spatial resolution of 127 km for this basin-wide ice 
model is incapable of resolving the narrow jet like features of the EGC. Also monthly 
ocean currents are used to drive the model which also reduces resolution in a temporal 
sense. 

C. ICE EDGE COMPARISONS 

The location of the ice edge boundary calculated by the model was compared with 
that determined by the Naval Polar Oceanography Center (NPOC). Sea ice analysis ice 
edge locations, identified on a weekly basis by NPOC, are based on a compilation of 
satellite imagery, aircraft reconaisance, and direct observations by reporting stations on 
land and sea (NPOC, 1986). The comparison that follows is restricted to the eastern 
Arctic to take advantage of its high climatic variability and the fact that it lay entirely 
within the ice model domain. The ice edge will be compared for 5 different time periods: 
11 March representing late winter; 20 May, the time of maximum ice extent; 8 July, the 
middle of the summer and time of maximum melt rate; 16 September, the time of mini¬ 
mum ice extent; and 25 November, the middle of autumn and a time of rapid ice growth. 

On 11 March (Figure 4.38) the ice edge extended from the east coast of Greenand 
at 65 °N to the western tip of Spitsbergen at 1S°N and then extended around the south¬ 
ern tip of the island at 75 °N, eastward to Novaya Zemlya, and then southward to the 
Kola peninsula. The modeled ice concentration for 11 March (Figure 4.39) appears to 
extend far south (an average of 180 km) of the actual ice edge line. The model ice edge 
appears to almost reach the northern coast of Scandinavia. The 50 percent concen¬ 
tration line appears to more closely match the real ice edge than any other ice concen- 
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tration level. One item of note is that the 40 percent and higher concentrations do 
recede back a little on the western side of Spitsbergen reflecting the presence of the warm 
influx of the West Spitsbergen Current indicating that the model is at least partially 
handling this heat flux correctly. The most likely explanation for the improper southerly 
projection of the ice edge is the use of climatology-derived oceanic heat fluxes which do 
not take into account the actual oceanic forcing. Also a factor is the large spatial re¬ 
solution (127 km) which can not resolve the ice edge which frequently is less than 50 km 
wide, especially during periods of on-ice winds. 

On 20 May the ice was at its greatest extent (Figure 4.40). The NPOC report indi¬ 
cated that the location of the ice edge was similar to that of 11 March, with the main 
difference being higher concentrations of ice inside the ice edge, especially between 
Spitsbergen and Novaya Zemlya. As in the previous case, the modeled ice edge (Figure 
4.41) appears to extend south of the actual ice edge by approximately 120 km east of 
Spitsbergen and 180 km west of Spitsbergen. The 90 percent concentration line seems 
to best approximate the real ice edge. 

By 8 July (Figure 4.42) the ice edge had begun to recede with the west coast of 
Spitsbergen and the southern coast of Novaya Zemlya becoming ice free. Open water 
was present above 80°A r near Spitsbergen and the ice concentration near the ice edge 
between Spitsbergen and Novaya Zemlya had decreased to 30 percent. The modeled ice 
edge (Figure 4.43) had also receded northward but was still south of the actual ice edge 
by 100 km. The model shows incorrectly that the west coasts of Spitsbergen and 
Novaya Zemlya are ice bound (at least partially) and the warm influx of the WSC is not 
observed. 

By 16 September the ice extent had reached its minimum. The ice edge depicted in 
Figure 4.44 starts at the eastern Greenland coast at 1A°N and extends northeastward 
above Spitsbergen over to Franz Josef Land and then eastward to the northern tip of 
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Novaya Zemlya. Spitsbergen and Novaya Zemlya are virtually ice free. Low concen¬ 
trations of ice extended deep into the Arctic. The model ice edge (Figure 4.45) was still 
70 km south of the actual ice edge but it compares better this time. Spitsbergen and 
Novaya Zemlya are correctly shown to be virtually ice free and ice concentrations are 
down throughout the Arctic basin. The model still has difficulty representing the warm 
intrusion of the WSC west of Spitsbergen. 

Figure 4.46 depicts the location of the ice edge on 25 November, examined for this 
period of rapid ice growth. The ice edge originated from the eastern Greenland coast 
at 65 °N and extended to the northwestern tip of Spitsbergen and ran northeastward to¬ 
ward the northern tip of Novaya Zemlya. It picked up again at the southern tip of the 
island and extended over to Siberia. For this period of ice growth the ice model (Figure 
4.47) correctly depicts the ice edge. It still tends to be too far south by 40 km but it is 
closer than the previous cases. The model seems to correctly reflect the warm heat flux 
of the WSC 

Overall, the ice model appears to do a better job of predicting the ice edge in the 
later part of the year (late summer,fall, early winter) than the first part of the year. 
Preller (personal communication) has indicated that PIPs behaves in a similar manner. 
Also of note is that PIPS is currently updated weekly with ice edge reports from NPOC 
and so never strays very far from the actual ice edge in the course of the week (Preller 
and Posey, 1989). 

D. DISCUSSION 

The main focus of this ice model study was to compare modeled ice drift to ground 
truth data using buoy trajectories in a qualatative manner. Ice drift and ice flow pat¬ 
terns were examined under a variety of environmental and forcing conditions at different 
locations. The entire ice field was not examined in detail; only that portion where buoys 
were present was studied (the exception being the Greenland Sea). Model performance 
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was studied under decreasing wind speed conditions, increasing wind speed conditions, 
steady state light or calm winds, steady high winds resulting from strong large pressure 
gradients, and gyre reversals. 

The model performed well in depicting the general ice flow patterns in the central 
Arctic basin. Here wind forcing exerted the greatest influence on ice drift. This was 
reflected in both the ice model and in the buoy reports. 

Some trends noted from the case studies are listed below. Under high wind speed 
forcing the ice model response was immediate and yielded an accurate depiction of the 
ice flow. Model speed differences varied 2 to 3 cm/s above and below the buoy drift 
speed. A 12 hour time lag (one time step) was noted for ice drift rates to slow during 
decreasing wind conditions. This probably resulted from an excessive build up of mo¬ 
mentum from the movement of ice. No time lag was observed for rapid increases in 
wind speed (resulting from a moving or developing storm). Model response was imme¬ 
diate under these conditions. However, running the model at a shorter time step might 
detect the presence of a time lag. The response to the summer reversal of the Beaufort 
Gyre was also immediate. 

Outside of areas of strong wind forcing the ice drift speed rapidly decreases, as in the 
19 June storm, where the only ice moving was that located near the low pressure center. 
Under light or calm wind conditions, generally under 10 m/s, the ice model drift is gen¬ 
erally very slow and highly variable in direction. Also at this level of forcing other fac¬ 
tors in the momentum balance equation exert more of an influence such as ocean 
currents and internal ice stress. Floes in areas of thinner and less concentrated ice 
tended to move faster than in areas of thicker and more concentrated ice. 

Some ice model differences with the buoy drift trajectories could be attributed to the 
inaccuracies in the NOGAPS pressure fields that drove the model. Instances were 
present where both pressure fields missed a small scale weather or oceanic circulation 
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feature indicated by the buoy drift pattern. Generally the pressure fields provided a 
reasonable depiction of the wind forcing. There were some trends in modeled ice drift 
direction as well. At low wind speeds the model had difficulty in assessing direction. 
Errors in excess of 45 degrees were noted. A possible reason is that oceanic forcing be¬ 
comes important at low wind forcing and the data base for this is monthly climatology 
with poor spatial resolution. At high wind speeds, on the other hand, the model per¬ 
formed well generally remaining within 10 degrees of buoy motion. Accuracy now de¬ 
pends on the accuracy (spatial primarily) of input pressure (i. e., winds) field. 

A comparison of model ice edge to NPOC ice edge reports was conducted for se¬ 
lected dates. It indicated that the model ice edge tended to be well south (> 100 km) 
of the actual ice edge through most of the year. Model performance was better in the 
latter part of the year than in the begining. Use of monthly oceanic heat fluxes with 
poor spatial resolution results in the loss of some detail in advected heat flux such as 
the West Spitsbergen Current. The 127 km resolution of the model grid will also result 
in some loss of detail along the ice edge. The PIPS model overcomes this by being up¬ 
dated by the weekly ice edge reports from NPOC, therfore maintaining a fairly accurate 
picture of the ice edge. Regional models along the ice edge currently being developed 
would also provide improved spatial resolution and depiction of the ice edge. 
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Figure 4.1 A time series oP monthly averaged sea ice extents Prom the available liSMR and 
SMMR data, 1973-1987, Por the entire Arctic polar region. The inset shows the average sea¬ 
sonal cycle over the years 1973-1976 and 1979-1986 (Parkinson and Cavalicri, 1979). 
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TIME STEP (12IIR INCREMENTS) 


l igurc 4.2 Model sea ice c.xi for the year 1986. Ice extent is any defined area 
that contains 15 percent or more ice. The x-axis is in time steps ol 12 hours for 
one year, and the y-axis is the ice extent in 10*£/n 3 . 
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figure 4.3 Model sea ice area for the year 1986. Ice area covers only the actual 
ice present within each cell. The x-axis is in time steps of 12 hours over a year, 
and the y-axis is the ice areal extent in 
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Figure 4.4 Model sea ice volume for the year 1986. lee volume is the product 
of the ice thickness and the ice concentration. The x-axis is in time steps o 1 ' 12 
hours Tor one year. The y-axis is ice volume in lO’Am’. 
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Figure 4.5 Buoy trajectories and buoy de¬ 
rived surface pressure for II March 1986. 
For estimation of buoy speed 1 cm equals 
23.1 cm/s. 
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Figure 4.G NOGAPS derived wind trajectories in m,s for 11 March I9RG The three vectors in the lower 
right corner represent speeds of 5, 10, and 20 m/s, respectively. 










Figure 4.7 Model ice drift trajectories in cm/s for 11 March 1986. The three vectors toward the bottom 
center represent speeds of 5, 10 and 20 cm/s, respectively. 



Figure 4.8 Buoy trajectories and buoy derived surface pressure for (a) 18 January, (b) 19 
January, (c) 20 January, and (d) 21 January 1986. For estimation of buoy speed I cm equals 

23.1 cm/s. 
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Figure 4.9 Model ice drift trajectories in cm/s for 20 January 1986. The three vectors toward the bot¬ 
tom center represent speeds of 5, 10 and 20 cm/s, respectively. 









Figure 4.10 NOGAPS derived wind trajectories in m/s for 20 January 1980. The three vectors in the 
lower right corner represent speeds of 5, 10, and 20 m/s, respectively. 
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Figure 4.11 Model ice drift trajectories in cm/s for 21 January 1986. 1'hc three vectors toward the 
bottom center represent speeds of 5, 10, and 20 cm/s, respectively. 












Figure 4.12 Buoy trajectories and buoy derived surface pressure for (a) 18 February, (b) 19 
February, (c) 20 February, and (d) 21 February 1986. For estimation or buoy speed I cm 
equals 2.11 cm/s. 
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Figure 4.13 Model ice drift trajectories in cm/s for 19 February 1986. The three vectors toward the 
bottom center represent speeds of 5, 10 and 20 cm/s, respectively. 








Figure 4.14 NOGAPS derived wind trajectories in m/s for 19 Februarv 1986. The three vectors in the 
lower right corner represent speeds of 5, 10, and 20 m/s, respectively. 









Figure 4.15 Buoy trajectories and buoy derived surface pressure for (a) 12 June, (b) 13 June, 
(c) 14 June, (d) 15 June, (c) l(i June, (f) June 1986. For estimation of buoy speed I cm equals 


23.1 cm/s. 


69 











Figure 4.16 Buoy trajectories and buoy derived surface pressure for (a) 26 July, (b) 27 July, 
(c) 28 July, and (d) 29 July 1986. For estimation of buoy speed I cm equals 2J.1 cm/s. 










Figure 4.16 Buoy trajectories and buoy derived surface pressure for (c) 30 July, (!•') 31 July, 
and (I) August 1986. For estimation of buoy speed 1 cm equals 23.1 cm/s. 
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center represent speeds of 5, 10 and 20 cin/s, respectively. 





Figure 4.18 Model ice drift trajectories in cm/s for 17 June 1986. The three vectors toward the bottom 
center represent speeds of 5, 10 and 20 cm/.s. respectively. 
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Figure 4.21 Buoy trajectories and buoy derived surface pressure for (a) 7 April, (h) 8 April, 
(c) 9 April, (d) 10 April, (c) 11 April, and (I) 12 April 1986. For estimation of buoy speed I 


cm equals 23.1 cm/s. 
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Figure 4.22 Model ice drift trajectories in cm/s for 7 April 1986. The three vectors toward the bottom 
center represent speeds of 5, 10 and 20 cm;s, respectively. 































Figure 4.25 Model ice drift trajectories in cm/s for 27 September 1986. The three vectors toward the 
bottom center represent speeds of 5, 10 and 20 cm/s, respectively. 
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Figure 4.26 Model ice drift trajectories in cm/s for 27 October 1986. The three vectors toward the 
bottom center represent speeds of 5, 10 and 20 cm/s, respectively. 



f igure 4.27 Buoy trajectories and buoy <1 
August, (c) 18 August, and (d) 19 August 
23.1 cm/s. 


ved surface pressure Tor (a) 16 August, (h) 17 
86. For estimation of buoy speed I cm equals 
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Figure 4.28 Model ice drift trajectories in cm/s for 16 August 1986. The three vectors toward the bot¬ 
tom center represent speeds of .5, 10 and 20 cm/s, respectively. 
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Figure 4.29 Model ice drift trajectories in cm/s for !9 August 1986. The three vectors toward the bot¬ 
tom center represent speeds of 5, 10 and 20 cm/s, respectively. 






I igurc 4.30 Buoy trajectories and buoy derived surface pressure for (a) 16 December, and (b) 
19 June 1986. l or estimation of buoy speed 1 cm equals 23.1 cm/s. 
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bottom center represent speeds of 5. 10 and 20 cm/s, respectively. 









Figure 4.32 Model ice drift trajectories in cm/s for 19 June 1986. The three vectors toward the bottom 
center represent speeds of 5. 10 and 20 cm is, respectively. 













Figure 4.33 Buoy trajectories aiul buoy derived surface pressure for (a) 20 May, and (b) 17 
September 1986. For estimation of buoy speed 1 cm equals 23.1 cm/s. 










Figure 4.34 Model ice drift trajectories in cm/s for 20 May 1986. The three vectors toward the bottom 
center represent speeds of 5, 10 and 20 cm/s, respectively. 











Figure 4.35 Model ice drift trajectories in cm/s for 17 September 1986. The three vectors toward the 
bottom center represent speeds of 5, 10 and 20 cmfs, respectively. 






Figure 4.36 Model ice drift trajectories in cm/s for 23 January 1986. flic three vectors toward the 
bottom center represent speeds of 5, 10 and 20 cm.s. respectively. 









Figure 4.37 Model ice drift trajectories in cm/s for l March 1986. The three vectors toward the bottom 
center represent speeds.of 5. 10 and 20 cm/s, respectively. 






















Figure 4.38 The NPOC eastern Arctic icc edge for 11 March 1986. 
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r igurc 4.40 The NPOC eastern Arctic ice edge for 20 May 1986. 
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.44 The NPOC eastern Arctic ice edge for 17 September 1986. 
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Figure 4.46 The NPOC eastern Arctic ice edge for 25 November 1986. 
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V. SUMMARY AND CONCLUSION 

A. SUMMARY 

As the U. S. Navy conducts operations in the Arctic, there will always continue to 
be a need for a way to predict the movement and behavior of the ice pack. Early com¬ 
puter models of ice prediction studied and forecast ice movement in isolation without 
considering the effects of the ice thickness, strength, and thermodynamic growth. Today 
the Navy uses the PIPS model, based on Hibler's (1979) dynamic-thermodynamic sea ice 
model, which combines a thermodynamic heat budget ice model with a dynamic ice 
model which is sensitive to the effects of ice thickness and ice strength interactions. 

The dynamic thermodynamic sea ice model simultaneously solves a system of four 
principle equations of momentum balance, ice rheology, ice thickness distribution, and 
ice strength. Momentum balance describes the change in ice movement or ice drift 
through time as the result of a summing of forces of wind and ocean stress, Coriolis 
force, a sea surface tilting term, and internal ice stress. Of these forces the wind stress 
term normally has the largest effect on the momentum balance. Ice rheology descibes 
the internal ice stress interaction through a viscous-plastic law that states that under 
light strain ice behaves like a highly viscous fluid, and under heavy strain it behaves like 
a plastic. The ice thickness distribution describes ice thickness and concentration 
through a seven layered scheme that parameterizes the effects of advection, divergence, 
growth and decay and includes a heat budget for the ice-ocean-atmosphere system. 
Finally, ice strength, which is used in the ice rheology, is described as a function of ice 
thickness and concentration. 

The ice model used in this study was adapted from one used by Owens and Lemke 
(1990) to model the ice cover of the Antarctic and Weddell Sea. For this study the ice 


103 







model was fitted with the PIPS Arctic grid (127 km resolution) for its domain and was 
driven by 1986 NOGAPS forcing data at 12 hour increments over one annual cycle. The 
forcing data included pressure fields for deriving geostrophic winds and atmospheric 
thermodynamic data to drive the heat budget. The monthly ocean currents and oceanic 
heat fluxes used by the ice model originated from the Hibler-Bryan coupled ice-ocean 
model. The ice model had an option for using a vertically varying prognostic oceanic 
mixed layer but for this study the mixed layer was held to a constant depth of 30 m. # 

To initialize the model it was started with an ice field of uniform thickness, concen- 
tration, and velocity and was driven by the 1986 forcing input data for three years to 
permit the ice field to stabilize. The third year results were examined and compared to 
buoy drift data to evaluate model performance under differing environmental conditions 

B. CONCLUSIONS AND RECOMMENDATIONS 

The main purpose of this thesis was to study and evaluate the performance of the 
NPS version of the PIPS Arctic sea ice model. The focus of the study was to evaluate 
the ice model's ability to depict the movement of ice in both speed and direction. Ice 
model results from 1986 were qualatatively analyzed and compared to Arctic buoy drift 
data under differing environmental conditions of decreasing wind speed, increasing wind 
speed, steady state light or calm winds, steady high winds resulting from strong pressure 
gradients, and gyre reversals. The case of the ice drift in the Greenland Sea was also 
examined. 

» 

It was found that the ice model generally performed well in depicting ice movement. 

Model performance was better at higher wind speed but performance at lower wind r 

speed was reasonable. Model reaction to changing weather and wind conditions was 
almost immediate. However, a 12 hour time lag was noted for decreasing ice drift 
speeds. There was less variability in ice drift direction with sustained winds above 10 
m/s Thinner and less concentrated ice tended to move at a greater rate than thicker and 
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more compact ice. Some discrepancies between model and buoy data could be traced 
to inaccuracies in the NOGAPS data used to drive the model. Model compari.ym with 
NPOC ice edge reports placed the modeled ice edge generally more than 100 km south 
of the actual ice edge, most likely resulting from the poor representation of ocean heat 
fluxes at the ice edge by the monthly oceanic heat fluxes. Weekly updates of the ice edge 
in PIPS as run at FNOC prevents this from being a problem. 

The sea ice model's reasonable depiction of ice motion warrants its continued use 
by the Navy as an operational model. Its chief limitation, correctly modeling the input 
data, will improve in time as more powerful computers come on line which will permit 
finer spatial resolution. The data sparse Arctic will always present a challenge with re¬ 
gards to small scale features contained within it. Regional models now under develop¬ 
ment will improve this aspect especially along the ice edge. 

Further areas of study would be to take the model output data and the buoy data 
and conduct a statistical analysis for a quantitative study. Additional studies could in¬ 
clude running the ice model with different forcing and parameters to see if there is an 
improvement in ice edge location and ice thickness accuracy, and running the model with 
different relaxation schemes to improve computational efficiency. 
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